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jil^' Abstract 

I The two-loop formfactor both for a U{1) x U{1) and a SU{2) x U{1) gauge theory with 

massive and massless gauge bosons respectively is evaluated at arbitrary momentum transfer 
q^. The asymptotic behaviour for ^ oo \s compared to a recent calculation of Sudakov 
logarithms. The result is an important ingredient for the calculation of radiative corrections 
to Z-boson production at hadron and lepton colliders. 



1 Introduction 

Precise measurements of cross sections for the production of massive and massless gauge bosons 
were one of the central topics of LEP experiments. At the LHC similar reactions, namely the 
production of W and Z-bosons, singly or in pairs, with or without additional quark or gluon 
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jets, will be crucial for precise studies of the electroweak and strong interactions. Single 
W- and Z-boson production will be used for the determination of parton distributions and 
eventually even for luminosity calibrations. A future linear collider, operating in the GIGA-Z 
mode, will measure the properties of the Z-resonance with unprecedented precision. All these 
measurements will rely on the theoretical knowledge of radiative corrections to better than one 
percent accurracy, perhaps even down to the level of several permille. QCD and electroweak 
radiative corrections, as well as their interplay, thus will be crucial for the interpretation of 
these results. 

QCD corrections to single W- and Z-production are identical to those for the Drell-Yan 
process and have been evaluated in two- loop approximation in [U [2] , those for Higgs boson 
production in |3j. Electroweak corrections for the on-shell process were computed some time 
ago (see e.g. [4j and references therein). The next step evidently requires to combine QCD 
and electroweak effects, resulting in non-factorizable terms of order OweakOs- For the inclusive 
Z decay rate these terms have been calculated for final states with up-, down-, and bottom- 
quarks O El [7] and turned out to be relevant for the precise determination of the strong 
coupling constant. However, these results cannot be directly applied to the production process 
and to more differential distributions. For the Z-boson such corrections for high px distribution 
have been obtained in [8j. 

In the present paper we describe conceptial developments and concrete results which are 
important ingredients for the complete evaluation of these non-factorizable terms of order 
«weakQ!s. In particular we consider those amplitudes which correspond to vertex diagrams 
with a virtual gluon attached to one-loop electroweak corrections. These are relevant for 
the "mixed" corrections of order Q;„eakQ:s to Z-boson production, and for hadronic Z decay. 
Essentially the same diagrams are also important ingredients for the combination of photonic 
and weak corrections to Z production in electron-positron collisions and to leptonic Z decays. 

Our study identifies the infrared singular as well as the finite parts, investigates the struc- 
ture of these singularities and shows how they can be combined with real radiation to arrive 
at a finite result. The infrared finite remainder will be presented in analytical form in terms 
of generalized poly logarithms. 

The form factor will also be investigated in the Sudakov limit M^/g^ <C 1. In the special 
case of an Abelian theory the result coincides with the one of [9] (see also |10j ) and allows to 
contrast the logarithmic approximation with the complete result. The calculational method 
relies on an approach that has already been successfully employed in a number of cases |lHll2j. 
General considerations restrict the structure of the final result to a sum of "basis functions" 
(in our case — generalized "harmonic" polylogarithms up to fourth degree) with specific argu- 
ments and prefactors. Calculating on one hand directly a large number of terms in the low 
expansion with the technique of large mass expansion, expanding the basis functions on the 
other hand, and equating the results, the coefficients in front of the basis functions can be 
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determined. In a final step most of the basis functions are transformed into Nielsen polylog- 
arithms, leading to a fairly compact result whose asymptotic behaviour can be analyzed in a 
straightforward manner. 

To facilitate the discussion, we present, in a first step, in section 2, the results for a U{1) x 
U (1) theory with one of the gauge bosons taken to be massive, the other one massless. The 
explicit analytical result confirms the factorization of the infrared singularities and allows to 
identify the infrared-finite remainder. In section 3 the formalism will be extended to a massive 
nonabelian theory and applied to the complete set of virtual corrections of order a^jeakOts, 
contributing to Z-boson production and decay. The triple-boson coupling leads to additional 
diagrams with additional generalized polylogarithms, which cannot easily be transformed into 
Nielson's polylogarithms. However, they can be evaluated numerically with high precision [13] 
and their asymptotic behaviour is under control. The paper concludes with a brief summary. 
Much of the formulae and calculational details will be collected in the Appendices. 



2 Abelian Theory 

For definiteness and simplicity we will, in a first step, consider the form factor in a ficticious 
U{l)xU (1) theory with one massive and one massless gauge boson and with coupling constants 
a and a' respectively. 

For the Abelian theory the form factor F will be defined as matrix element of an external 
current 

7^F{q,M) = I dxe-'^m^'lMxM. (1) 

Here ip and ip' denote on-shell massless fermions of momenta p and p' = p + q, respectively, M 
the mass of the gauge boson. 
In a perturbative expansion 

F{a,a',q,M,e)= ^ f£j f|-j (g, M, e) (2) 

m,n=0 \ / \ / 

one needs to evaluate the expansion coefficients /(™'"-). In Born and one-loop approximation 
they are given by 

= 1, (3) 

^(10) 7 2 2 + 32, , ,2(1 + 2)2/ vr^N 

/^'•"^ = --- - + ^^\ogi-z)+ {Li2il + z) - —) , (4) 

where z = q^ /M"^ + iO, Cn = C(^) is the Riemann (^-function and the infrared singularities are 
controlled by dimensional regularization in d = 4 — 2e dimensions. In the euclidean region 
q^ < 0, so that no imaginary parts appear in the above formulae. 
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The two-loop result for the massless case, f^^'^\ can be found e.g. in |14tll5j. The two-loop 
result for the fully massive case, f^'^'^\ is only known in the large limit [16]. The evaluation 
of the mixed corrections is drastically simplified by the fact that the infrared singularities 
factorize within infrared evolution equation approach [17^ [TB] , which gives in our case 

F{a,a' /M^e) = F^^,,i^ss(.a' ,q,e)F{a,a' ,qyM^) , (6) 

with -Fmassless = J^i'^' f^''^)"' f^^'^^Q^ ^) denoting the formfactor for the massless theory and F 
being free from infrared singularities. The function F can again be expressed as double series, 
and the coefficients depend on the ratio z = q^ /M"^ only. The terms = fi'^fi) coincide by 

definition with those valid for the massive U (l)-theory. The evaluation of the nonfactorizable 
part of the two-loop contribution 

F(1'1) = 4>{q^/M^) (7) 

will be the central result of this section. 

The Feynman diagrams necessary for this computation have two thresholds: at = 
and at q^ = M"^. The analytical structure of vertex diagrams of this type has been explored 
in [TT]. The coefficients of an expansion in g^/M^ (and M'^ /q^) can always be expressed as 
combinations of so-called harmonic sums [18] or more generally — nested harmonic sums [19j . 
These sums correspond to (generalized) polylogarithms (pQj) [21] of arguments itg^/M^ and 
their generalizations — harmonic polylogarithms [22] (see also |23]). This structure suggests 
the following method for the evaluation of Feynman integrals. First, using the method of large 
mass expansion |24) . one calculates a large number of coefficients of the series in q^/M'^. From 
the basis functions (polylogarithms) one then constructs an Ansatz with unknown coefficients 
Xi- Equating Ansatz and series one obtains a unique answer for parameters Xj. This method 
has been applied earlier |251lll| to various scalar vertex masterintegrals. (In a different context 
the method has also been applied in [26j). Here it is applied to amplitudes deduced from a a 
set of realistic Feynman diagrams representing a physical process and leading to amplitudes 
with irreducible numerators and shrunken lines. 

A few comments on this procedure are in order. First, the main problem is to write down 
the correct prefactors in the Ansatz. Empirically one finds that the presence of a numerator 
or the absence of a line may lead to the additional factors M"^ /q^ or (M^/g^)^ in front of 
polylogarithm^. Therefore such factors should also be included in the Ansatz. Second, only 
five functions could not be represented as Nielsen polylogarithms with the argument q'^ /M"^. 
These remaining functions belong to the class of harmonic polylogarithms [23 discussed in 
more detail in the Appendix. 

^ In a series representation such multiplications lead to shifts of the summation index in c„. Indeed, if z = /M"^ 
then, e.g. j Y.n=i Cn^" = ci + Y.n=i Cn+i-z" and so on. 
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Instead of expanding the amplitude in /M"^ one could find the differential equation (see 
|27j ) for a diagram and again apply an Ansatz based on polylogarithms. This approach has 
recently been used for similar two- loop vertex diagrams in [28]. 

Altogether 16 one-particle-irreducible two-loop vertex diagrams contribute to the formfac- 
tor. These diagrams can be obtained from the one-loop one shown in Fig. [1^ by adding one 
gluon line. The two-loop, one-particle reducible diagrams which are obviously products of 
one-loop diagrams contribute to the term aa' f^^'^^ f^^''^^ and are not repeated here. We also 
do not display the contributions to the fermionic wave function renormalization, which receives 
contributions from additional 6 diagrams. For the generation of the input the program DIANA 
[29j has been used, for the evaluation and expansion a program written in FORM [30]. The 
evaluation of the Dirac traces has lead to about 700 different integrals. For most of them the 
asymptotic expansion was performed up to order 45 which required in total several hours of 
CPU time on a Pentium IV processor. For the remaining, most complicated cases (nonplanar 
diagram) up to 60 expansion coefficients had to be computed. For this purpose the parallel 
version of FORM [31j, running on an SGI machine with multiprocessor SMP architecture, was 
used. 

The function </>(z) can be cast into the following form (here and below z = /M"^ + fO) 
- ^ \ ' (6L2 + 24(2 - 24C3) log(l + z) + (-4L2 _ 6L - 2OC2) log^(l + z) 



z^ 

+^ log3(l + z)L + 81og2(l + z)U2{-z) - 121og(l + z)U2{-z) 
-161og(l + z)U^{-z) + 161og(l + z)Si,2(-2) - 16Li2(-z)C2 - AU2{-z)L^ 
-8Li^(-2) + 16Li3(-z)L - 24Li4(-z) - 12Si,2(-^) + 16Si,2(-z)L + 16Si,3(-z) 

-16S2,2(-^) + 24/i(z) - 48i/3(z) + m2{z) + 32Hi{z)j 

\ ^ 2>z + z"^ ( 

+ -32Li2(z)C2 - 8Li2(z)L2 - 8Li2(z)Li2(z) + 32Li3(z)L 

-48Li4(z) + 32S2,2(-z)^ + (^721og(l - z)C2 + 181og(l - z)l^ 

+361og(l - z)Li2(z) + 36Li2(z)L + 72Si,2(z)^ + ^-^(32(2 + I2C3) 

34 + 51z^ 16 + 23Z o 2(1 - z)(13 + 27z) ^ ,^ , 4(3 + 4z)^., , 

L H V - '-\ log(l - 2;) + ^ — 5 — -\.\2{-z)L 

z z z^ 

2(l + z)(3 + 5z) n^^T^j-f ^ . 4(11 + 9z) ^. , , 
^ -{log{l + z)L + Li2{-z)] + ^ -Li2{z) 

4(3 + 22-3^2) 4(9 + 42_6z2) 2(8 - z) 
-2 ^Li3(-z) - ^ ^2 ^Li3(z) + ^— ^, (8) 

where L = log(— g^/M^), = ((a) is the Riemann (-function, S'a,fe(z) are Nielsen polyloga- 
rithms [20j . The functions /i(z) and Hj{z) are defined and discussed in Appendix A. 
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3 Z-Production 



For definitenes and simplicity, consider, in the next step, Z-boson production in quark-antiquark 
annihilation. To fix the notation, we recapitulate the one-loop results. The weak corrections to 
the Born term can be split into those involving the exchange of W- and Z-bosons, (FigH^a)) 
and those involving the triple-boson coupling (Fig(T](b)). The combination of photonic and 
QCD corrections follows essentially from the two-loop QED or QCD results and will not be 
addressed here. 

For a light quark the form factor can be decomposed as follows 



At the Born level the expressions for the form factors Fji and Fj^ are given by 

6 6 

Fr = i — gn , Fl = i — gi , 



(9) 



(10) 



with g^i = —Qs']y/cw and gL = {h — Qs]y)/cw being the right- and left- handed couplings 
of a quark to the Z-boson. Here I3 is the third component of the isospin of a quark, Q its 
electric charge and sw = sin 9\y and cw = cos 6w denote sine and cosine of the weak mixing 
angle, respectively. 

Including radiative corrections and adopting a form similar to eqlHthe formfactors can be 
cast into the following form 



Fn = ^— + 

su/ V 47r 



sw 

e 

I — 

sw 

a 



gR + 



gL 



^T^slr 



9RPA{q^/m\, m\) + CF 



9R(pA{q'^/ml 



V3 I 2 I 2 2 \ , 9L / 2 / 2 2 \ , ^3 / 2 / 2 2 \ 

9LPA{q imz^mz) + —pA{q /m^^rny^) + cw—p^A{q /m^r,mw) 



+ Cf 



a, a 



(11) 



where first factors in the brackets in the above equations represent the QCD corrections. Terms 
given by pA and pNA account for the one-loop electroweak corrections. The abelian part pA 
is defined by the diagram of the abelian type (Fig. la) and obviously closely related to f^^''^^ 
defined in eqlH The unrenormalized resuHo is given by 



Pa{z,M'^ 



ln(MV/i^ 



+ 



2 + 3z 



z 



z 



log(-z) + 



2(1 + z)^ 



z^ 



vr 



Li2(l + z)--^j. (12) 



^We shall not discuss issues related to renormalization, since the non-factorizable part, which is the quantity of 
interest in this paper, is independent of the renormalization scheme. 
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The nonabelian part /9na(^5 M2) receives corrections from both diagrams of Fig. la and Fig. lb. 
It is given by 



A{z,M^) 



3 

e 



2z'' z 



with 



^ = ln 



v/1 - 4/z + 1 + iO 



(13) 

(14) 

(15) 



y/l-4/z -1 + iOj 

The function A(z) can be taken from [33] (see also [3ll[35l|36] and references therein for one- loop 
calculations in the Standard Model). We do not include the terms from the renormalization 
of the coupling and the Z-boson wave functiorlfl which follow from textbook prescriptions and 
will not be considered in this work. 

Evaluated for arbitrary q"^ 7^ M|, the above results are gauge dependent and are presented 
in Feynman gauge. For the offshell case they can be considered as building blocks for a complete 
calculation. 

The functions </>a(-z) and (/'na(-2), representing the non-factorizable terms of 0{aaf), are 
written in a form completely analogous to the electroweak one-loop terms. The function <Pa{z) 
has been given in the previous section. The nonabelian part (ptssAiz) involves new functions — 
generalized poly logarithms. Our result in Feynman gauge reads 




(b) V 



Figure 1: Diagrams, contributing to the vertex Zqq (a) and (b). The two-loop diagrams are obtained 
by attaching one virtual gluon in all possible ways. The case (b) represents nonabelian part. That 
gives contribution 0na(-2) in the text. Diagram (a) with W exchange also contributes to 0na(2;). 



<Pa{z) = (16) 
</'na(2) = -24>a{z) 

+4^^ + 48^^i/o,-r,-r,-l(-^) - l2H.r,-r,-l{-z) + 8//o,-r,-r (-^) 



■^Hence the function A(z) considered in differs by subtracting the term 3/e — 3 ln(Af ^//i^) — 1/2. Furthermore, 
a typo in [33] has been corrected, flipping the signs of the terms proportional to I and 



7 



+4 -2 H_r,-r{-z) - ^ '-gH^r,-l{-z) 

-2^ ^2 '-gR^r{-z) + ^ (8C2 + 2L2j - 4 Li3(z) 

66 + 49Z, ri - z)(13 + 34z) , , , 5 + 9z ^ 

+2 Lis z - 4^ ^ log 1 - z - 16 L 

z z-^ z 

+ ^ ^ ( 961og(l - z)C2 + 241og(l - 

+48 log(l -z)Li2(z) +48Li2(z)L + 9681,2(2) j 
'^^^^^2^^'^ ( -32Li2(z)C2 - 8Li2(z)L2 - 8Li2(z) 



z2 



(/)a(1 + iO) = 14 + 72C2/2 - 64C2/i - —l\ + 22C2 - 28C3 + I6C4 - 128Li4(i: 



+32Li3(z)L - 48Li4(z) + 3282,2(2) j • (17) 

The function (/>na(-2) receives contributions not only from digrams of Fig. 1(b) but also from 
those of Fig. 1(a) with the exchange of 1^-boson. The functions i?,.. are considered in more 
detail in Appendix B. For the special case = one finds 

16, 

Y 

32 

+i7r(85 + 32/2 + 24/2 _ q ^ 14^2 _ ^^20(3) (18) 

o 

= -2.1073 - 19.0331i, (19) 

32 70 1 84 

(/>na(1 + iO) = -16 - 144(2^2 + 128C2^i + y /| + yC2 + — Cs - 236(4 

+26^ + 256Li4(i) - 84-^Ls2(f ) - y^Ls2(f ) + 96(Ls2(f ; 

64 

+m(54 - 64/2 - 48/i + — /^ - 28C2 + 48C3) (20) 
= -7.5880 + 16.7194i , (21) 

with I2 = log 2. 8ubstituting the actual masses of the W- and Z-bosons {z = m|/m^ = 
1.2856) we find: 

(/>a(1.2856 + iO) = -1.3598 - 30.4095i , (22) 

(/)na(1-2856 + iO) = -10.1248 + 35.0336i . (23) 

In the limit q'^ ^ 00 the function (p, as given by Eq. ([8|) coincides with the result of [9] 
where the power unsupressed logarithmic and constant part have been evaluated. For the 
leading and the first power suppressed term we find 

</)(z) = (3-24C2 + 48C3)log(-z)-2 + 40C2-84C3 + 14C4 

+ - f (-26 + 8C2) log'(-2) + (-120 - I6C2 + I28C3) log(-z) 



z 



-I88-8C2-8C3 + II6C4) +0 ( ^ ) (24) 
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In Fig. [2] the exact result is compared with the Sudakov approximation and with the approx- 
imation including the first power-suppressed term. For electroweak interactions the mass of 
the gauge boson can be taken to be of order 100 GeV, the characteristic energy of order one to 
two TeV. For one TeV the relative error of the Sudakov approximation (the logarithmic plus 
constant term) amounts to 15%, at 2TeV it is reduced to 2.5%. 



1 ' 1 ' 1 ' 
~f(z) 


1 ' 1 ' 


1 


exact 

Sudakov approximation 
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- 
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\\ 
\ 

\ 


\ 

\ 
\ 


I.I.I, 


I.I. 


1 



-400 -300 -200 -100 



Z 



Figure 2: Non-factorizable two- loop correction to the abelian formfactor in the euclidean regime 
[z = g^/m^). The solid line represents the exact result, the dashed line the Sudakov approximation 
and the dash-dotted line includes the power suppressed terms. 



4 Conclusions 

Using the technique of asymptotic expansions and the knowledge of the general structure of 
integrals we evaluated analytically the two-loop formfactor in a U{1) x U{1) theory with one 
massive and one massless gauge boson. In the Sudakov limit full agreement is abtained with 
[9], where the logarithmic and constant terms had been evaluated obtained. We furthermore 
perform the same caldulation for a SU{2) x U{1) theory and derive the non-factorizable part 
of the two-loop formfactor in the Standard Model. As an application we evaluate the mixed 
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virtual 0{aas) radiative correction for Drell-Yan production of the Z-boson. 
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5 Appendix A 



In this Appendix we consider the asymptotic behaviour of the most complicated basis functions 
in the limit z = q^/m? oo. Most of basis functions can be expressed in terms of Nielsen 
polylogarithms and then the standard transformations formulae can be applied to go from 
argument 2; to 1/z (see |20^ 121]). Therefore we will consider here only the five special cases, 
mentioned previously, where complications arise. 

In our calculation the following functions appear in addition to usual Nielsen polyloga- 
rithms: 

h{z) =i?_i,o,i(^) , 
Hi{z)=H^i^o,iAz), 

H2{z) =i?_i,o,0,l(2;) , 
/73(^)= ^-1,-1,0,1 

i74(2;) =i7o -1,0,1(2:) , 

where Ha^ . 4{z) are harmonic polylogarithms defined in [22J . 

These functions correspond to the alternating Taylor series in z: 

Kz)=-f:'-^^^{-zT, 

1 77- 
n=l 

H,{z) = -± ^-'''^''-'\ -z)\ 

n=l 

00 



n=l 



00 



^ 5_3(n - 1) + S'_2,i(ra - 1) - 5'i(n - l)5_2(n - 1) 

^3(2) = - 2^ {-z) 

n=l 



00 



n=l 



with finite harmonic sums Sa{n) = Z]j=i 1/j" and S-ain) = Z]j=i(~l)''/j" a-iid S^2,i{'n) = 
X]j=i(— l)-'5'i(j)/j^. It is interesting to note that the function Hi cancels in the final result 
([H]) for the formfactor but is present in the particular integrals. 
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Following [TT] it is not difficult to write down simple integral representations for the above 
series, e.g. 

z 

Kz)= f-^U,{x), (1) 

J 1+ X 


z 

^i(^)= / TT^Si'2(x), (2) 



z 

^2(^)= / if^Lislx), (3) 



dx 



H^{z)=log{l + z)h{z)- J -^Li2(x)log(l + x), (4) 



2 

H,{z) = log{z)h{z) - / Li2(x) log(:E) . (5) 



Now the integrals can be expressed in terms of Nielsen polylogarithms of nonlinear ar- 
guments and only one harmonic polylogarithm function H2 (this choice being not unique, 
however). We have 

Hz) = ^Si,2(^') - Si,2(^) - Si,2(-2) + ln(l + z)Li2(z), (6) 
Hiiz) = log(l + 2)Si,2(^) + ^Si,3(z=^) - Si,3(-z) + l^z), (7) 



Hsiz) = log(l + Z) (iSi,2(z2) - Si,2(z) - Si,2(- 



where 



+ i log2(l + z)U2{z) + \s,,s{z') - Si,3(z) - h{z), (8) 

Hi{z) = ^S2,2(^') - S2,2(^) - S2,2(-^) + log(l + z)Us{z) - H2{z), (9) 



$(z) = -^C4 + ^log^slog^+ ^log2s(Li2(s) -Li2(-s) 



- log s(Li3(s) - Li3(-s)) + Li4(s) - Li4(-s), 

with s = (1 - z)/{l + z). 

In order to find the asymptotic behaviour for z ^ oo one needs to use the standard formulae 
for polylogarithms and for the function H2{z) = //_i^o,o,i(-2) the inversion formula (A. 6) from 
[32j . It is important to take care of imaginary parts, therefore we approach the cut in g^-plane 
from above, which means that z is replaced by z + iO. Thus we obtain 

h{z + iQ) = -^log3z + 2C2logz-^C3 + ^(^-^log^^-logz + 2C2^ 

+ .vr{ilog^.-ic2 + ilog. + ^}+o(l), (10) 

11 



13 57 

Hliz + iO) = — log^Z- -C2l0g^^ + C3l0g2 + — C4 



+ i Q log3 z + i log2 z - 3(2 log ^ + Ca - 3C2 - l) 

+ iTT log^ Z + C2 log 2 - ^CS + ^ ( - ^ log^ Z - logZ + C2) I + 0(^)(,11) 

F2(z + i0) = -^log^Z + C2log2^-^C4 

+ 1(^-1 log^ 2 - \ log2 z + (2C2 - 1) log Z + 2C2- 2j 

+ ,vr|^log3z-^C3 + ^Qlog^^ + logz + l)} + o(^), (12) 

Hsiz + iO) = -^log^Z + C2log2^-^C3logZ-^C4 

+ i(^-ilog3z + (2C2 + l)logz-^C3 + l) 

+ ^vr{ilog3.-lc2log. + ^C3 + ^Qlog^.-^C2-l)} + o(^), (13) 

Hiiz + iO) = -^log^^ + C2log^^-^C3log^ + ^C4 

+ ^Qlog2z + 21ogz-2C2 + 2^ 

+ zvr|^log3z-lc2logz + ^C3 + ^(-logz-2)|+o(^). (14) 

Finally we used the program hplog [13] to check numerically the asymptotic behaviour of the 
iJ-functions. 



6 Appendix B 

In this appendix we consider the //-functions contributing to the nonabelian part of the form- 
factor. For the definitions and recursive constructions of these functions we refer to |28j . 
However, for completeness we give here explicitly the definitions of the functions which appear 
in our calculation. The following six new functions arise in the evaluation of the two-loop 
nonabelian formfactor: 

Z t2 

dt2 f dti 





Z t2 

dto f dt 



/ X /■ f dtl 



12 



r,— r, — 1 



dU 



t3 



dU 



t2 



y/tsits + i) J y/t2it2+i) J 1 + ^1 



dti 



dU 



t2 



dU 



0,— r,— r, 



-l(^) 



17 



t4 



dU 



dh 



t2 



dti 



(18) 
(19) 
(20) 





In the formula (|17p for the nonabehan part the functions with odd number of indices "— r" 
appear always with the factor 



1 



(21) 



It is easy to check that H^r and H^^-i cannot be expanded in the Taylor series of small 
arguments (they have a branche point at zero), but the combinations gH-r and gH^r-i can. 

The integral representations given above are not very suitable for the analysis and numerics. 
The ultimate task would be to relate them to the usual (harmonic) polilogarithms. In order 
to do this one should choose a "right" variable. From the previous expirience ^JLj it is known 
that for the diagrams, posessing a branch point at = 4m^, the appropriate variable is given 
by {z = /m?) 



l + ^z/{z-A) 
In terms of y the y-factor (j2ip is expressed as 

1-y 



9{-z) 



1 + y ' 



and the required //-fucntions take form 



H, 



0,—r,—r 



H-r-r-l{ — Z 
HQ-r,-r-l{ — Z 



= - log y , 
= 2 log y, 

= \ log^ y + ^Li2(-y3) - Li2(-y) - ^C2 , 

= -]:\og^ y + log(l - y) log^ y - 2Li3(y) + 2 log y Li2(y) + 2C3 , 
o 

= -\ log3 y + logy + ^Cs + Li3(-y) - ^Li3(-y3) , 
b 6 6 9 

11 2 89 1 

= ^ log"^ y - gC2 log^ y - gCa log y + — C4 - Li4(-y) + — Li4(-y^) 

+ 2Si,3(l - y) - |c2Li2(l - y) + 21n(l - y) f ^Cs + Li3(-y) - ^Li3(-y 



(22) 
(23) 

(24) 
(25) 

(26) 

(27) 

(28) 
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where 



+ 2Li2(y)( Li2(-y) - ^Li2(-y3) ) + 2iVi(l) - 2iVi(y) , 



7Vi(y) = j'^U2{t)Hl-t + t% iVi(l) = -ll^, 



(29) 



(30) 
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As it is seen from the above fomulae, the i?-functions with index "— r" can be rewritten in 
terms of harmonic polylogarithms but of nonlinear argument y. 
In the hmit when z — ^ +00 + iO we obtain 



(31) 



g{-z)H_ri-z) = log ^ + i (21og z - 2) + m |-1 - ^1 + C? Q 

gi-z)H^r,-i{-z) = \ V - y C2 + (log' z - 2 log z - y C2 - 1 
+i7r |- log z - ^ (21og z - 2)| + O (^^ 



+ 117 { — log Z + 



z 



i?-.,-.(-.) = ^log'.-3C2-^ 

Ho,-r,-r{-^) = \ log^ ^ - 3C2 log Z + 2(3 + ^ (2 log Z + 2) 



+ 



+Z7r<{ -^log'2; + C2 



+ 



1 



z" 



H-r,-r,-l{—z) 



- log'' Z - — C2 log Z+-C3 + - 
6 3 3 z 



,2 20^ , 
log2 Z + —C2 + I 



+iTr 



1 



log^ z + -C2 + - log z 



z^ 



15 2 

i?0 -r -r = ^ log^ 2 - 3C2 log' Z + -(3 \og Z + 7(4 



+^ log' z + (2 + ^C2 - ^Cs) log 2; + 1 



20 



1 o 4 2 

- - log 2:+ -C2logZ- -C3 

2 log z - 2 - |C2 + ICs - 2 log zH + O 



+i7r 
1 

"^zV ' ~ 3'" ' 3" 

And finaly we give the values of i7-functions at the particular point z = \: 



9{-m 

g{-m- 

H. 

Ho-r. 

H—r - 



0,—r,—r. 



-1) 
-1) 

-1) 



3^/3' 
2 Ls2(f) 
3 V3 ' 



yC2, 

-l) = tC3-|^Ls2(f), 



-1) 
■1) 



:C3, 



■^C4+|(LS2(f))\ 



where („ is the Riemann (-function and Ls„(a;) is the log-sine integral defined as 

X 

Lsn{x) = - J log"-i(2sin^)dt. 



(32) 
(33) 

(34) 
(35) 



(36) 

(37) 
(38) 
(39) 
(40) 
(41) 
(42) 



(43) 
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In particular the constant Ls2(|), sometimes denoted as Clausen's integral Cl2(|) (see, e.g., 
[2Tj). is given by 

Ls2(f ) = 1.014941606409653625 . . . (44) 
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